function Fig2_c_1()
data = 48;
Fcolor='k';
Bcolor='w';
Fsize=18;  %font size

result = CS1_0_extraction_2D_Phi(data);

assignin('base','dataOneHalf',[result.az_r;result.az_phi])

%%%% plot azimuthal average %%%%
g=figure('name',['Grid:',num2str(data)]);

phi_line = 1000* (result.az_phi - result.avg );
errorbar(result.az_r,phi_line,1000*result.az_phi_err./(result.phi_grouped_count').^0.5,'LineWidth',1);
% plot(result.az_r,phi_line);

g.Color = Bcolor;
ax = gca;
ax.Color=Bcolor;
ax.XColor=Fcolor;
ax.YColor=Fcolor;
xl=xlabel('r (nm)');
yl=ylabel('(meV)');
ax.FontSize=Fsize;
yl.FontSize=Fsize;
xl.FontSize=Fsize;
axis([0,125,-4.2,1])

end